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The idea that the parton system created in relativistic heavy-ion colhsions (i) 
emerges in a state with transverse momenta close to thermodynamic equihbrium 
^ ' and (ii) its evolution at early times is dominated by the 2-dimensional (transverse) 

, hydrodynamics of the ideal fluid is investigated. It is argued that this mechanism 

I may help to solve the problem of early equilibration. 
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1. It is now commonly accepted that evolution of the partonic system created 
I in heavy-ion collisions at RHIC energies is best described by hydrodynamics of 

an almost ideal fluid [1]. In particular, particle transverse-momentum spectra 
and asymmetry of the transverse flow, as represented by the parameter V2 [2], 
are reasonably well reproduced by the hydrodynamic approach. 

At the same time it is also realized that the hydrodynamic picture encoun- 
ters a rather serious challenge: it requires very early thermalization of the 
system. This follows from the fact that the asymmetry of the transverse flow 
is produced most effectively at the very early stage of the evolution (when 
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the pressure gradients are largest). To obtain the asymmetry parameter V2 
consistent with data it is necessary to start the hydrodynamic evolution at 
the time below 1 fm after the collision takes place. But application of hydro- 
dynamics demands, as the necessary condition, the local equilibration of the 
system. Such fast equilibration is not easy to achieve with elastic perturbative 
cross-sections. This puzzle of early thermalization was also widely discussed 
and several exotic mechanisms were proposed for its solution [3.4] . It seems 
fair to say, however, that no one was yet accepted as fully satisfactory. 

In the present paper we explore the possibility that, at its early stages, the 
hydrodynamic evolution applies only to transverse degrees of freedom of the 
partonic system created in high-energy collisions. The main reason of interest 
in such a study is the observation that the early equilibration, if at all possible, 
is particularly difficult to achieve in longitudinal direction. This is because an 
elastic collision does not change significantly the direction of the colliding par- 
tons and thus it requires very many interactions to produce a locally isotropic 
distribution from the initially strongly anisotropic one 0- 

At the same time one may argue that the equilibration of the partonic spec- 
trum in transverse direction can be obtained much easier. Indeed, it is well 
known that the transverse momentum spectra observed in nucleon-nucleon col- 
lisions are well described by the Boltzmann distribution . As these spectra 
reflect - at least to some extent - the distribution of produced partons [7f8|9] . 
it is clear that the partonic system produced in hadronic collisions emerges 
already in a state close enough to equilibrium in the transverse direction so 
that no much more action is needed to achieve the goal, see e.g. [TU|lllll2j . 

At this point it is important to emphasize that such a description may be 
adequate only during a certain time after the collision. Indeed, as the time 
goes on, the interpartonic interactions will tend to equilibrate the system also 
in longitudinal direction and thus the standard, 3-dimensional hydrodynamics 
will probably take over. 

To study the possibility of such purely transverse hydrodynamic evolution and 
to have also a connection with physics of the relativistic heavy-ion collisions 
we discuss a simple model implementing these ideas. The parton distribution 
is constructed as a superposition of many clusters whose longitudinal motion 
satisfies the Bjorken in-out condition t] = y where 



is the spatial rapidity of the cluster and y is its rapidity. Transverse momenta 
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in each cluster follow the equilibrium distribution in the local rest frame of 
the fluid. 



This is realized by the following Ansatz for the parton distribution function 

F{x,p) = /|| g{r,r];x^,p_i),lans (2) 
where f\\ describes the longitudinal motion of clusters, 



S(y — T]) 

/|l = noS{p^t - Ez) = no , If par (3) 

Tm± 

with T being the longitudinal proper time, —t^ — z'^, uq is the normalization 
factor giving the density of clusters per unit of rapidity, and g represents 
the distribution of transversally equilibrated partons. In our study we neglect 
parton masses and take for g the Boltzmann distribution 



p U^^ 

= exp ( — ^ ) , Idisg (4) 



T 

where is the four velocity of the fluid, 

— {uQCOsh.y,Ux.iUy,UQSmh.y) — (woCOshTy, WosinhTy), (5) 

with u'^ being the 4-velocity of the fluid in the frame where the longitudinal 
momentum vanishes 



u'^ — (i^o, Uy, 0) = {uq, u±, 0), ul — u]^ — l.lumu (6) 



Using these ideas, the hydrodynamic equations were derived from the Ansatz 
(??) - (??) and the transverse expansion of the fluid was studied numerically 
for various initial and final conditions. 

Starting from the initial profile determined by the density of the wounded 
nucleons inside the colliding nuclei, we have found that 

(i) With the proper choice of the initial temperature Tj it is possible to ob- 
tain parton spectra that are consistent with the experimentally observed pion 
transverse- momentum distributions. 

(ii) The calculated elliptic flow parameter V2 agrees with data, provided the 
flnal temperature Tf (i.e. the temperature of transition from 2-dimensional 
(2D) to 3-dimensional (3D) regime) is taken substantially higher than the 
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expected freeze-out temperature. This gives a rather short time available for 
the 2D evolution. 

(iii) This short evolution time implies that the transverse size of the system 
is relatively small, smaller than that obtained from the HBT measurements, 
thus consistently leaving space for further, 3D expansion of the system. 

We thus conclude that the 2D hydrodynamics may indeed be a reasonable 
description of the partonic system during the first few fermis after the collision 
and it gives a good starting point for a realistic treatment. If confirmed by more 
detailed studies, this would imply serious changes in the present understanding 
of the quark-gluon plasma evolution. 

These conclusions differ substantially from those reached in a study of a similar 
problem by Heinz and Wong [13] . Using a different implementation of the ideas 
advocated in the present paper, they found that it is not possible to obtain 
a satisfactory description of the elliptic flow and concluded that the partonic 
system must be necessarily in a 3D equilibrium from the very beginning of 
the evolution. Our example shows that this is not necessarily the case. In the 
last section we discuss the possible reasons for this discrepancy. 

In the next section we discuss the general form of the energy-momentum tensor 
following from the Ansatz (??) - (??) and the corresponding hydrodynamical 
equations. The numerical analysis of the evolution of the system are presented 
in Section 3. The last section summarizes our conclusions. 



2. To obtain the hydrodynamic equations we first derive, using the Ansatz 
(??) - (??), the general form of the particle density and entropy density 4- 
vectors and of the energy-momentum tensor. Somewhat lenghty but straight- 
forward calculations give (with h = c = 1) 

dy—^p^'fw g = no Ug- — t/^, Idensn (7) 
zvr r /vrr 



where Vg is the number of internal degrees of freedom {vg = 16 for gluons) and 



V = {smhri,0,0,coshr]) (9) 

is the 4-vector defining the longitudinal (z) direction. It is also interesting to 
evaluate the entropy flow, giving 



4 



= -no Ug / rfy— ^p'^/ii 9{\ogg - 1) = rio Ug^U^.ldenss (10) 
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From these results we also read that the energy, entropy and particle number 
densities in the local rest frame are 

dE dN 1 dS T2 

'^o'^q — ; -Ti r = r = T^o^qi:, — • (11) 



d'^x±dz TTT d'^x±dz 3 d'^x±dz 2tit 

where the equation of state dE/{d'^x±dz) = 2P was used. Note that depen- 
dence of these densities on temperature is weaker than that known from 3D 
thermodynamics. 

Using dz = rdr] we thus obtain the densities per unit of spatial rapidity 



dE _ T\ _ dN _ _1 dS _1 

d'^x±dr] ^ IT ' ^ d'^x±dT] ° ^2tt 3d'^x±dri 3 ^' 

(12) 



which are identical to the densities proper for 2D thermodynamics. 

The hydrodynamic equations are obtained from the energy-momentum conser- 
vation laws df^T'^'^ = which also imply the conservation of entropy d^S^ = 0. 
Entropy conservation implies 



dr[uoT^] + V± ■ [n±T^] = 0, legs (13) 
while the momentum conservation provides two other equations 



dr[u±T^uo] + [T^ui]iV± ■ Ui_) + (m± ■ Wi_)[T^ua] + V±TV3 = 0./egt(14) 

These are three equations for the three unknowns: the temperature T, and 
two independent components of the 4- velocity and Uy. Their characteris- 
tic feature is that they do not depend explicitly on the space-time variables 
{T,x,y). Thus, they do not select any special point in space time. We have 
solved these equations numerically with the method described in [14| which is 
a direct extension of that proposed in [Lb]. The details of the procedure are 
published in [T6j . 

3. To study the evolution of the system following from (??) and (??) we need to 
specify the initial conditions, adequate for the physical situation encountered 
in Au—Au collisions at RHIC energy. To this end we assume that the profile (in 
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Fig. 1. Transverse- momentum spectra of positive pions measured by the PHENIX 
Collaboration in the centrality class 30-40% (solid line) |T8] and the model spectra 
of gluons for various choices of Tj and Tf. The lowest dash dotted curve corresponds 
to Ti = 200 MeV and Tf = 180 MeV. The highest dash dotted curve was obtained 
for Ti = 300 MeV and Tf = 180 MeV. The four almost parallel lines represent our 
results for Tj = 250 MeV and for four final values of the temperature: Tf = 200 
(long dashed line), 180 (dash dotted line), 160 (dashed line) and 140 (dotted line) 
MeV. Note that the original tt"*" experimental spectra collected at ^snn = 200 GeV 
have been multiplied by a factor of 3 to account for the total hadron multiplicity. 

transverse coordinates) of the initial energy density €2 is given by the density 
of participants. This density is determined, for a given centrality, from the 
Glauber formulae. The following discussion is given for centrality 20 — 40%, 
corresponding to impact parameter b ^ 7.9 fm. The initial temperature Tj at 
the origin, x± = 0, is taken as a free parameter. 

The parton spectra are evaluated using the Cooper-Frye prescription [17] 



where E is the surface at which the 2D evolution comes to the end. In our 
exercise we have taken S to be the surface of constant temperature. 

The parton transverse-momentum spectrum evaluated from the model de- 
pends substantially on the initial temperature of the system but is practically 
insensitive to the final temperature (this observation was already made in 
|13j). In Fig. 1, the spectrum obtained for various values of initial tempera- 
ture Tj is compared with measured spectrum of pions [IB]. One sees that the 




(15) 
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correct slope [fj is obtained for Tj ^ 250 MeV. This initial temperature is used 
in all subsequent calculations. The normalization, which is arbitrary in our 
model, was calculated with tiq = 1. 

At the fixed initial temperature, results for elliptic fiow are sensitive to the 
value of the final temperature Tf. This is shown in Fig. 2 where V2 calculated 
from the model is plotted versus p± for various final temperatures. One sees 
that, as expected, V2 grows with decreasing Tf (at a fixed Tj, smaller Tf means 
longer time of evolution and thus more time to develop the fiow). As already 
observed in pS] , the absence of longitudinal pressure implies that the obtained 
values of f 2 are higher than those expected from 3D evolution (with the same 
initial conditions). To bring V2 close to the observed one [19], a fairly high 
Tf 180 MeV is needed. This agrees with the point of view formulated in the 
first section: The final temperature should not be interpreted as the "freeze- 
out" temperature but rather as the temperature when the 2D character of 
the equilibrium changes into the 3D one. This may happen well before the 
freeze-out and hadronization. So the high value of Tf is not surprising. 

The fact the f 2 grows with time indicates that our approach can be adequate 
only at the initial stage of the evolution of the system. After this initial stage 
other effects such as the early transition to 3D hydro with a reduced sound 
velocity or with viscous effects, hadronization, and switching to transport 
description will reduce the growth of V2 shown in Fig. 2. 

Having fixed the initial and final temperatures, we can also estimate the life- 
time of the 2D evolution. It turns out to be of about 4 fermi (for the 20 — 40% 
centrality class). 

4. In conclusion, we have investigated the consequences of the hypothesis that 
the partonic system produced in a collision of two heavy ions is created in 
a state close to thermodynamical equilibrium in transverse direction, while 
its longitudinal structure is characterized by a freely-streaming collection of 
clusters. Two observables were studied using the Cooper-Frye formula at the 
surface of constant temperature with the initial energy density profile deter- 
mined from the density of the number of participants at a given centrality. 

(a) The slope of the transverse momentum spectrum was adjusted to the one 
measured for the pion spectra [TS] giving the initial temperature at the center 
of the system Ti = 250 MeV. 

(b) The calculated elliptic fiow parameter V2 gives the final temperature (at 



^ At this stage the experimental shape of the spectrum is not correctly reproduced 
since we do not include the effects such as, e.g., resonance decays and hard scatter- 
ing. 
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Fig. 2. The elliptic flow coefficient V2 as the function of the transverse momentum. 
The PHENIX experimental results for pions and kaons in the centrality class 20-40% 
and for the collision energy ^/snn = 200 GeV [19j are compared to the model cal- 
culations with Tj = 250 MeV and for four final values of the temperature Tj = 200, 
180, 160 and 140 MeV. The curves are denoted in the same way as in Fig. 1. The 
best agreement is obtained for Tf ~ 180 MeV. 

which the transition from 2D to 3D regime takes place) Tf = 180 MeV, 
substantially higher than the expected hadronization temperature. The time 
needed for the 2D evolution was estimated to be much shorter than that 
normally needed to achieve the freeze-out. This confirms the point of view 
that the equilibrium changes from 2D to 3D well before the freeze-out and 
hadronization. 

We thus conclude that the initial period of 2D transverse equilibration and 
hydrodynamic evolution of the parton system is helping to solve the problem 
of early equilibration, most pertinent difficulty of the present application of 
hydrodynamics to physics of heavy-ion collisions. It would be certainly inter- 
esting to investigate this possibility in more detail. 

Several comments are in order. 

(i) Our conclusions differ from those obtained in [13]. As far as we can see, 
apart from certain technical details, there are essentially three reasons for this 
discrepancy. 

First, the density distributions implementing the assumption of free- streaming 
in the longitudinal direction differ by a factor depending on the transverse 
mass. This leads to rather serious consequences. Our ansatz (2) - (4) implies 
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that at fixed rapidity the evolution of system obeys the rules of truly 2D ther- 
modynamics of an ideal fluid, whereas in [13] the temperature dependence of 
the energy and number densities follows that of 3D thermodynamics and the 
resulting hydrodynamic equations should be interpreted as an effective de- 
scription of the viscous hydrodynamics. Our energy-momentum tensor is more 
symmetric and, consequently, the hydrodynamic equations differ from those 
used in [13]. We stress that our approach is thermodynamically consistent 
and the conservation laws for the energy and momentum, d^T'^'^ = 0, lead to 
the entropy conservation, dfj^S^ = 0, in the way very much similar as in the 
standard relativistic hydrodynamics. 

Second, in Ref. [13] a new timescale parameter tq is introduced. This parameter 
- absent in our formulation - plays an important role in the argument used 
in [13] to reject a two-dimensional evolution. The large values of tq, found in 
the fitting procedure by Heinz and Wong, are in their opinion inconsistent 
with the very idea of transverse thermalization because such thermalization 
may take place only at the beginning of the evolution of the system, when 
the considered times r are small (with tq < r). The replacement tq — no/Tn± 
(which formally transforms their model into our approach) solves this problem 
if no is interpreted as the density of 2D clusters in rapidity. 

Third, the authors of [13] confront their results for the elliptic flow with the 
earlier 3D hydrodynamic calculations, whereas we fix our parameters by com- 
parison with the present data. Our results shown in Fig. 2 indicate that we 
can reproduce the experimental values of f 2 already at high values of the tem- 
perature where most of the effect is created. The 3D hydro results shown in 
Fig. 7 of [13] are above the data, hence they cannot be used as the reference 
point for rejecting the concept of transverse thermalization. 

(ii) In our calculations we have assumed a sharp transition between the 2D and 
a possible 3D evolution. This simplification, necessary to obtain a reasonably 
tractable problem, can be removed in future, more sophisticated analyses. This 
raises interesting questions about the nature of this transition and about the 
final fate of the clusters. 

(iii) It should be clear that our investigation represents only the first step 
towards a fully realistic description. Many details, as e.g. the mechanism of 
2D — > 3D transition, the relative distribution of clusters, their internal (longi- 
tudinal) structure, the rapidity dependence of the system, are left for future 
work. Nevertheless, we feel that we produced some compelling arguments for 
the existence of a substantial period of purely transverse thermodynamic equi- 
librium and hydrodynamic evolution of the partonic system (in the form of 
the ideal fluid) at the very beginning of the collision process. 

Acknowledgements 



9 



Discussions with W. Czyz, B. Mueller, St. Mrowczyiiski, and K. Zalewski 
are highly appreciated. This investigation was partly supported by the MEiN 
research grant 1 P03B 045 29 (2005-2008) and Polish Ministry Of Science and 
Higher Education grants Nos. N202 153 32/4247 and N202 034 32/0918. 



References 

[1] For a review see, e.g., E.Shuryak, Nucl. Phys. A750 (2005) 64. 

[2] J. Y. Ollitrault, Phys. Rev. D46 (1992) 229. 

[3] St. Mrowczyiiski, Acta Phys. Pol. B37 (2006) 427. 



[4] M. Strickland, talk at Quark Matter 2006, |hep-ph/0701238 
[5] F. Becattini, Z. Phys. C69 (1996) 485. 



[6] F. Becattini and U. Heinz, Z. Phys. C76 (1997) 269, |hep-ph/ 9702274| 



[7] A. Bialas, M. Bleszynski, and W. Czyz, Nucl. Phys. Bill (1976) 461. 
[8] A. Bialas and W. Czyz, Acta. Phys. Pol. B36 (2005) 905. 
[9] A. Bialas and A. Bzdak, Phys. Lett. B649 (2007) 263. 



[10] A. Bialas, Phys. Lett. B466 (1999) 301, |hep-ph/9909417 



[11] W. Florkowski, Acta Phys. Pol. B35 (2004) 799, nucl-th/0309049 



[12] S. Steinke and J. Rafelski, J. Phys. G32 (2006) S455, |nucl-th/0607066[ 



[13] U. Heinz and S.M.H. Wong, Phys. Rev. C66 (2002) 014907, and U. Heinz, 
private communication. 



[14] M. Chojnacki, W. Florkowski, Phys. Rev. C74 (2006) 034905, |nucl-th/0603065 



[15] G. Baym, B. L. Friman, J. P. Blaizot, M. Soyeur, W. Czyz, Nucl. Phys. A407 
(1983) 541. 

[16] M. Chojnacki and W. Florkowski, Acta Phys. Pol. B39 (2008) 721, 
rarXiv:0710.5871 [nucl-th]. 

[17] F. Cooper and G. Frye, Phys. Rev. DIO (1974) 186. 

[18] S. Adler et al., PHENIX Cohaboration, Phys. Rev. C69 (2004) 034909, 



nucl-ex/0307022 



[19] S. Adler et al., PHENIX Collaboration, Phys. Rev. Lett. 91 (2003) 182301, 



nucl-ex/0305013 



10 



